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Abstract The effects of initial tree spacing on wood density 
at breast height were examined for 22-year-old Japanese 
larch ( Larix kaempferi). The experiment involved the use 
of three plots with different initial tree spacing densities 
(300, 500, and 1000 trees/ha). For five trees selected from 
each plot, the total tree height, diameter at breast height, 
height to the base of the live crown, and crown diameter 
were measured. Ring width and wood density for individual 
growth rings were determined by X-ray densitometry. A 
mixed-effects model was applied for fitting the radial varia¬ 
tion in wood density in relation to initial spacing. Models 
having various mean and covariance structures were tested 
in devising an appropriate wood density model. The model, 
consisting of the mean structure with quadratic age effects 
and heterogeneous first-order autoregressive covariance, 
was able to describe the radial variation in wood density. 
Closer spacing of trees (1000 trees/ha) resulted in a faster 
increase in wood density from the pith outward than for 
more widely spaced trees, indicating that initial tree spacing 
may influence the age of transition from juvenile to mature 
wood. 

Key words Japanese larch • Silvicultural treatment • Initial 
spacing ■ Linear mixed model • Wood density 


Introduction 

Japanese larch ( Larix kaempferi ) shows rapid growth and 
good adaptability to cold climates, and has been planted in 
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subalpine regions and the northern area of Japan for several 
decades. 1 The area afforested with L. kaempferi accounts 
for about 30% of all plantation forest areas in Hokkaido, 
northern Japan. 2 Wood products from the immature stands 
were low in quality owing to the large volume of juvenile 
wood, and their primary use was limited to piles, civil engi¬ 
neering lumber, pallets, and packing materials. However, 
as the stands matured, the timber has been recommended 
for use as construction material, including structural glued 
laminated timber. Consequently, Japanese larch has good 
potential as an important wood product for sustainable 
harvest from intensively managed plantations. 3 ' 4 

Numerous studies have examined the effects of silvicul¬ 
ture practices on wood properties. 5 Of these, tree spacing 
(encompassing initial planting density and thinning treat¬ 
ment) is one of the most common silviculture practices for 
controlling the growth and form of individual trees, and has 
been studied intensively. Wood density was found to 
decrease following thinning treatment of Pseudotsuga men- 
ziesii, 6 Pinus banksiana, 7 Pinus radiataf and Picea abiesf 
but in other studies, thinning had minor or no effects on the 
wood density of Pinus taeda , w12 Pinus ponderosa, 13 and 
Picea rubens. 14 It has also been reported that thinning 
increased wood density in Pinus resinosa , 15 Larix occiden¬ 
talism and Pseudotsuga menziesii - 17 In Japanese larch, Koga 
et al. 18 found that heavy thinning significantly affected the 
growth rate and tracheid length at breast height, but had 
little effect on wood density. There are some reports con¬ 
cerning the effects of the initial spacing on wood properties 
in Pinus taeda and Pinus ellioti, 19 Picea ahiesf' Pinus resin- 
osaf 1 and Pinus banksiana ; 22 however, results are not 
consistent. These inconsistencies could be due to species 
variation, genetic variation, climatic factors, precision of 
measurement, and sampling techniques. 5-23 ’ 24 

Wide initial spacing of trees and thinning treatments 
decrease competition among trees, improving light avail¬ 
ability and access to water and nutrients in the soil, thus 
increasing plant growth. 5-25 Relationships between growth 
rate and wood properties have also been investigated in 
many studies, but the results have been controversial. 
Growth rates affect wood properties at different plant ages. 
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and ignoring this effect could result in growth rates being 
incorrectly identified as the cause of differing wood proper¬ 
ties. 26,27 In black spruce (Picea mariana ), Zhang 28 found that 
the negative correlation between ring density and ring width 
tended to be weaker with increasing age, and suggested that 
fast growth has a less negative effect on wood density with 
increasing tree age. Similar trends have been reported in 
Abies balsamea 29 and Larix species. 30,31 Thus, it is essential 
to consider initial tree spacing, and the optimum tree age 
and time for thinning treatments in the management of tree 
plantations during a rotation age. 5,23 

Sequences of yearly measurements of wood characteris¬ 
tics are considered as longitudinal (time series) data, 32 and 
the analysis of longitudinal data is able to evaluate the age 
effect on variation of these characteristics. In this case, 
these data are repeated measures made of the same 
characteristic on the same observational unit (tree, disk, 
and ring), and such data generally present temporal auto¬ 
correlation, heteroscedasticity, and nonstationarity of the 
mean. 24,32,33 Objections are frequently raised to statistical 
analyses that ignore correlation among observations due to 
bias in the variance of estimated coefficients and hypothesis 
tests. 34 Gregoire et al. 35 introduced the mixed-effects analy¬ 
sis technique as a means of accounting for correlation 
among observations. Correlations among observations 
made on the same subject or experimental unit can be 
modeled using random effects, random regression coeffi¬ 
cients, and through the specification of a covariance struc¬ 
ture. 36 Recent studies have developed models predicting 
variation in wood properties within a stem or cross section 
using the mixed-effects analysis for many species. 37 43 In this 
article, we develop a model describing the radial variation 
in wood density in relation to initial spacing and discuss the 
effect of spacing on wood density based on the model pre¬ 
dicted value. 


Materials and methods 

Experimental site and sampling 

The trees used in this study were obtained from an experi¬ 
mental site (43°14' N, 143°33' E; approximately 220 m ele¬ 


vation) in the Ashoro Research Forest of Kyushu University 
(Hokkaido, northern Japan). The annual average precipita¬ 
tion and temperature in the research forest are 782 mm and 
5.9°C, respectively. Japanese larch seedlings were planted 
in 1984 at an initial plant density of 2000 trees/ha. When the 
stand was 9 years old (1992; stand density 1000 trees/ha), 
three plots (each 625 m 2 ; 25 x 25 m) were established in the 
stand. Two treatment plots were thinned (plot A: 300 trees/ 
ha; plot B: 500 trees/ha), and the third plot (plot C) was 
retained as a nonthinned control. At the time of treatment, 
crown closure had not occurred in the stand, and conse¬ 
quently the effects of the thinning treatment were similar 
to the effects of initial spacing. The stand density at 9 years 
old was relatively low, but the reason for this was unknown 
(perhaps vole browsing). 

In June 2005, when the stand was 22 years old, five trees 
were randomly selected from each of the three plots. The 
diameter at breast height (DBH) of each selected tree was 
measured, and the tree was harvested. Total height (TH), 
height to the base of the live crown (HB), and crown 
diameter (CD) were measured (Table 1). A 5-cm-thick, 
knot-free sample disk was excised at breast height (1.3 m) 
from each harvested tree, and the disks were transported to 
the Hokkaido Forest Products Research Institute for 
processing. 

Wood density measurements 

A 2-mm-thick (longitudinal direction) strip that included 
the pith was cut from each sample disk along the arbitrary 
radial direction. Resin was extracted from the strips for 1 
week in a solution of benzene-ethanol (2:1), and the strips 
were then dried to an equilibrium moisture content of about 
12 %. 

Wood density and ring width measurements were 
obtained using X-ray densitometry. 4445 Although this tech¬ 
nique can provide various wood density and width charac¬ 
teristics for individual growth rings (e.g., earlywood density, 
latewood density, and latewood proportion), the mean 
wood density within ring and annual ring width data were 
used in this study. Each strip was scanned from the pith 
toward the bark in two radial directions, and the wood 
density and width for individual rings were expressed as the 


Table 1. Characteristics of the sample trees at harvest 


Plot 

NTA (trees/ha) 


TH (m) 

HB (m) 

DBH (cm) 

H/D ratio 

CD (m) 

A 

300 

Mean 

18.75 

6.41 

27.55 

68.23 

3.46 



Min 

17.24 

4.80 

25.47 

65.38 

2.73 



Max 

20.27 

8.30 

30.78 

74.22 

4.25 



SD 

0.98 

1.18 

2.07 

3.16 

0.49 

B 

500 

Mean 

18.60 

8.06 

25.66 

73.42 

3.05 



Min 

18.05 

6.12 

22.92 

57.86 

2.65 



Max 

19.19 

9.42 

31.20 

83.34 

3.61 



SD 

0.49 

1.08 

2.87 

8.54 

0.33 

C 

1000 

Mean 

18.90 

11.06 

20.33 

95.01 

1.96 



Min 

16.50 

10.30 

15.60 

76.88 

1.63 



Max 

20.30 

11.90 

25.89 

107.54 

2.48 



SD 

1.49 

0.57 

3.79 

11.13 

0.31 


NTA, Number of trees per unit area after treatment; TH, total tree height; HB, height to base of live crown; DBH, diameter at breast height; 
H/D ratio, tree height divided by diameter at breast height; CD, crown diameter; SD. standard deviation 
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mean of both values. In order to determine whether the 
wood density pattern varied with ring position in terms of 
cambial age, wood density of individual rings were averaged 
into 5-year groups. The number of annual rings included in 
the sample disks was 19 or 20, and thus the total number of 
observations was 60 (15 trees x 4 groups). 


Statistical analysis 

Observations from the same unit (tree, disk, and growth 
ring) tend to be correlated. As the wood density measure¬ 
ments were made ring by ring at breast height, it is reason¬ 
able to assume that they would be temporally correlated. 24 33 
The mixed-effects analysis technique is frequently used to 
account for correlations among observations. 35 46 47 In the 
mixed model, the correlation structure of the dependent 
variable can be taken into account by allowing the param¬ 
eters to vary randomly from one individual to another 
around the fixed population mean. 48 A linear mixed model 
for expressing the wood density of a specific tree is: 

\Y i =X i p + Z i b i + e i 

Ibj ~ N ( 0 , D) (1) 

[£, ~ MO, X,) 

where Y [ is the n, dimensional vector of density observa¬ 
tions for tree i, 1 < i < N, N is the number of trees, X t and 
Z, are (n, x p) and (n, x q) dimensional matrices known 
covariates,/I is the p dimensional vector containing the fixed 
effects, bj is the q dimensional vector containing the random 
effects, and e, is a n, dimensional vector of residual com¬ 
ponents. D is a general (q x q) covariance matrix with (/, j) 
element d» = d„ and Z, is a (n l x «,) covariance matrix 
that depends on i only through its dimension n h which 
means the set of unknown parameters in Z, will not depend 
upon L 

The model has two types of error terms: the b i are the 
random subject effects and the e, are observation level error 
terms. In the linear mixed model we assume only that the 
b t and e, are independently distributed with zero mean and 
variances D and Z, respectively. On the basis of this model, 
the observations vector T, has a multivariate normal distri¬ 
bution with an expected value of 

E (Y t ) = Xffi (2) 

and variance 


In Eq. 4, a is the vector of all variance and covariance 
parameters found in V h which means that a consists of the 
q(q + l)/2 different elements in D and of all parameters in 
Xi- If ce was known, the maximum likelihood (ML) estima¬ 
tor p, obtained from maximizing Eq. 4, conditional on a , is 
then given by: 47 




(5) 


and its variance-covariance matrix then equals: 
var (j= [|^var(lj)»'A: i ] 

( 6 ) 

= ( 7 ) 


where W t equals V) 4 ( a ). If a is not known but an estimate 
a is available, it is possible to set V t = Vj(a) = IT) 1 , and 
to estimate p by using the Eq. 5, in which IT, is replaced 
by TT ; . 

A common method for estimating a is the restricted 
maximum likelihood estimation (REML). This method 
enables estimation of the variance components, explicitly 
taking into account the loss of the degrees of freedom 
involved in estimating the fixed effects. The REML estima¬ 
tor for the variance component a is obtained from maxi¬ 
mizing the likelihood function of a set of error contrasts 
U: 


U = K'Y 


Y\~Y 


Yn-2-Y_ 

Uv-i-eJ 


( 8 ) 


where K is an \n, x (n [ - p)\ full-rank matrix with columns 
orthogonal to the columns of the X matrix. The vector U 
then follows N((),K l V(a)K). The likelihood function of the 
error contrasts can be written as follows: 49 


N V 2 1 V - 1 / 2 (V 

L{a) = {2K)- (n - p)j2 J j X]X i x] rxtvr'x, ni^r 1/2 

1=1 i =1 i=l 

xw{~f i (Yi-X,P)'vr 1 (Y l -X l p) 


(9) 


var {Y i ) = V i =Z,DZ] + Z i 


(3) where 


where superscript notation t indicates the transposed 
matrix. 

The classical approach to inference is based on estima¬ 
tors obtained from maximizing the marginal likelihood 
function (Eq. 4) with respect to 0 , which is the vector of all 
parameters in the marginal model for Tj. 36 


L ML (0) = n{(2^r /2 |Y i (a)r V2 x 

1=1 L 

exp(-i(I) - Xjpy Vr\a){Yj - *,/»))} 


(4) 


a = —f J (Y i -X i P) l V i r\a)(Y-X i p) (10) 

n-i-p tS 

It should be noted that closed-form solutions for the REML 
estimators of variance components can be obtained only in 
particular cases. 36 In most cases the restricted log-likelihood 
must be maximized in an iterative way such as the Newton- 
Raphson algorithm. 50 

The statistical analyses and model parameter estimations 
were performed using the MIXED procedure of the SAS 
V8 software (SAS Institute Inc.). 51 Under the linear mixed 
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model (Eq. 1), the data vector Y t for the z'th tree is assumed 
to be normally distributed with mean vector Xfi and covari¬ 
ance matrix Vi = ZDZ\ + Therefore, fitting of linear 
mixed models implies that an appropriate mean structure 
and a covariance structure need to be specified. Selection 
of the appropriate mean and covariance structures were 
based on maximum-likelihood criteria, Akaike’s informa¬ 
tion criterion (AIC), 52 Bayesian information criterion 
(BIC), 53 and the Likelihood Ratio Test. The AIC and BIC 
were calculated as: 

AIC = -21n(L)+ 2*7 (11) 

BIC = -2 ln( L) + <7 ln( A*) (12) 

where L equals the maximized ML or REML likelihood 
function, q is the number of variance components in the 
model (i.e., the number of parameters in a), and N* is N 
(total number of observations) in case of ML estimation 
and N - p in case of REML estimation (p is the number of 
fixed parameters). The AIC and BIC were used to discrimi¬ 


nate among models having the same fixed effects (vari¬ 
ables), but different covariance structures. 

The likelihood ratio statistic (-2 In A v ) was defined as: 36 


-21nA v = -21n 




(13) 


where L ml equals the likelihood function (Eq. 4), and 
0MLO and 0 VII , are the maximum likelihood estimates under 
the null hypothesis and the alternative hypothesis, 
respectively. 


Results and discussion 

Radial trends in growth and wood density 

The radial (cambial age) trends from pith to bark for ring 
width and wood density measurements at breast height in 
the trees from each spacing group are shown in Fig. 1. The 


Fig. la, b. Radial trends in 
a ring width and b wood density 
at breast height for each spacing 
plot. Thick solid lines , mean 
value of each plot. Symbols 
represent each five individuals 
within the plot 
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trends in ring width for plot A (300 trees/ha) and plot B 
(500 trees/ha) were quite similar (Fig. la); the values 
increased with cambial age up to about 10 years of age and 
then declined gradually. The ring width in the control trees 
(plot C; 1000 trees/ha) showed the same trend as those for 
the other two plots up to about 10 years of age, but then 
decreased markedly with further age. 

Although the effects of initial tree spacing on the radial 
growth rate were quite evident, the effects on the wood 
density were unclear. Wood density for the control trees 
(plot C) decreased to 10 years of age and then steadily 
increased with further age (Fig. lb). In contrast, wood 
density for trees in plots A and B were almost stable with 
age and showed no clear increase in radial variation. 


Model development 

The longitudinal (time series) data on wood density were 
analyzed using linear mixed models (Eq. 1) to determine 
the effects of initial spacing on the radial trends. In model 
development, an appropriate mean structure was first 
selected from the next three models. The first preliminary 
model (model 1) is the saturated mean structure, incorpo¬ 
rating a separate parameter for the mean response at each 
5-year group within each spacing group: 

Y„=Po + (Pi SA; + ftSB,. + ftSC;) 

+ (ftSA,. + ftSB; + p 6 SC,) x AGE,y (=5) 

+ (ft SA, + ft SB, + ftSC,) x AGE,y (=10) (14) 

+ (ft 0 SA, + ftiSB; + ft 2 SC,) x AGE,y (=15) 

+ (PuSAj + ft 4 SB, + ft 5 SC,) x AGE 1 y (=20 ) 

where Y, is the density observation of the 2 th tree, measured 
at y'th age (AGE,y). SA„ SB,, and SC, are indicator variables 
defined as 1 if the ith tree was obtained from Plot A, Plot 
B, and Plot C, respectively, and 0 otherwise. P are the 
parameters for each fixed effect. 

Second, the model was simplified by assuming a linear 
trend between the observation and age (model 2): 

y„ = ft + (ftSA, + ftSB, + ftSC,) 

+ (ftSA, + ftSB, + ftSC,) x AGE,y 1 ; 

From the radial trends shown in Fig. lb, third model is 
assumed to have a linear as well as a quadratic age effect 
within each spacing treatment (model 3): 

Y,j = ft + (ftSA,. + ftSB,. + ftSC,) 

+ (ftSA, + ft SB, + ft SC,) x AGE, (16) 

+ (ftSA,. + ftSB,. + ft SC,) x AGE* 

AGE was expressed in decades rather than years as in the 
original dataset (i.e. 0.5, 1.0, 1.5, and 2.0 years). This was 
done to avoid the possibility that the random slopes for the 
linear and quadratic time effects show too little variability, 
which might lead to divergence of the numerical maximiza¬ 
tion routine. 36 

The results of fitting three different mean structures to 
the wood density data are shown in Table 2. It was assumed 
that these three models have an unstructured covariance 


Table 2. Summary results for fitting three different mean structures to 


wood density 
structure 

data for 

models with 

an unstructured 

covariance 


Par 

~2L ml 

-21nAjv 

df 

p 

Model 1 

16 

-315.0 




Model 2 

7 

-294.9 

20.1 

9 

0.017 

Model 3 

10 

-306.4 

8.6 

6 

0.197 

Model 3-2 

7 

-305.4 

9.6 

9 

0.384 


Par, Number of parameters; minus twice the maximized ML 

log-likelihood value; -2 In X N , likelihood ratio statistic testing vs model 
1; df, degree of freedom; P, P-value corresponding to the likelihood 
ratio test statistic 


structure. The likelihood ratio statistic (-21 n A v ) is the dif¬ 
ference between -2 L ml of the current and reference model. 
The ML likelihood ratio test comparing model 2 with model 
1 rejected the null hypothesis of linearity (P = 0.017), and 
thus the saturated mean structure should not be simplified. 
In contrast, the quadratic model (model 3) was well fitted 
and was acceptable compared with model 1 (P = 0.197). 
Furthermore, the fixed effect of spacing was not significant, 
based on the type III tests 51 (null hypothesis H 0 : ft = ft = 
ft = 0; P = 0.569; data not shown). The quadratic model 
without spacing effect (model 3-2) was also acceptable 
compared with model 1 (P = 0.384) and model 3 (P = 0.801, 
data not shown). Hence, the mean structure of model 3-2 
(Eq. 17) was used for further inference and simplifying of 
the covariance structure. 

y„=Po+( ftSA,. + ft SB, + ft SC,-) x AGE, 

+ (ftSA; + ftSB,. + ft SC,) x AGE* 

The profile of the ordinary least squares (OLS) residuals 
(r, = y, - Xfi ols ) and the variance function of these residuals 
are shown in Fig. 2. This OLS estimator p OLS is consistent 
for p from the theory of generalized estimating equations 
(GEE). 54 The OLS residuals varied greatly for each tree, 
and the variance was not stable over time (heterogeneous). 
The estimated covariance matrix obtained from fitting of 
model 3-2 is: 


f 0.0018 0.0010 0.0010 0.00073 
0.0010 0.0010 0.0010 0.0008 
0.0010 0.0010 0.0011 0.0008 
,0.0007 0.0008 0.0008 0.0009, 

and the corresponding correlation matrix is: 

fl.OOO 0.733 0.742 0.5933 
0.733 1.000 0.953 0.853 
0.742 0.953 1.000 0.845 
,0.593 0.853 0.845 1.000, 

There was no clear evidence for a particular trend in the 
matrices, although the covariance between two measure¬ 
ments showed a slight decrease as a function of the time lag 
between them. Hence, the following six models were fitted 
to explore the covariance structure; (1) unstructured, (2) 
random intercepts, AGE and AGE 2 , (3) random intercepts 
and AGE, (4) random intercepts, (5) heterogeneous 
Toeplitz, denoted by TOEPH, and (6) heterogeneous first- 
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Table 3. Summary of results for fitting of various covariance models to wood density data for 
models with a quadratic mean structure without spacing effects (model 3-2) 


Model 

2^reml 

AIC 

BIC 

df 

P 

Unstructured covariance 

-249.7 

-229.7 

-224.5 

10 


Random intercepts, AGE and AGE 2 

-241.0 

-227.0 

-224.5 

7 

0.034 

Random intercepts and AGE 

-235.3 

-227.3 

-226.5 

4 

0.025 

Random intercepts 

-230.3 

-226.3 

-226.0 

2 

0.013 

TOEPH 

-240.8 

-226.8 

-224.3 

7 

0.031 

ARH (1) 

-240.5 

-230.5 

-229.2 

5 

0.101 


-2LJJEML, Minus twice the maximized REML log-likelihood value; AIC, Akaike's information 
criterion; BIC, Bayesian information criterion; df, number of parameters (degrees of freedom) 
in the corresponding covariance model; P, P-value corresponding to the likelihood ratio test 
statistic compared with the unstructured covariance structure; AGE, tree age; TOEPH, hetero¬ 
geneous Toeplitz; ARH (1), heterogeneous first-order autoregressive 




Fig. 2a, b. Profiles of the ordinary least squares (OLS) residuals (a) 
and the variance in OLS residuals (b) for wood density at breast height. 
Open triangles , plot A (300 trees/ha); open circles , plot B (500 trees/ha); 
filled squares , plot C (1000 trees/ha) 

order autoregressive, denoted by ARH (1). The last two 
covariance structures assume the correlations to be a func¬ 
tion only of the time lag between the measurements. 36 In 
addition, the Toeplitz structure assumes no relation among 
all correlations in different bands. The ARH (1) structure 
assumes all correlations to be of the form p", where p is the 
correlation coefficient between two measurements and u is 
the time lag between the measurements. Table 3 shows the 
results of fitting of these covariance models to wood density 
data for models with the quadratic mean structure (model 


Table 4. Restricted maximum likelihood estimates and standard errors 
(SE) for all fixed effects and covariance components in model 4 


Effect 

Parameter 

Estimate 

SE 

Intercept 

ft 

0.4437 

0.0185 

AGE x spacing effect 

SA 

ft 

-0.0671 

0.0270 

SB 

ft 

-0.0327 

0.0270 

SC 

ft 

-0.0968 

0.0270 

AGE 2 x spacing effect 

SA 

ft 

0.0221 

0.0104 

SB 

ft 

0.0090 

0.0104 

SC 

ft 

0.0402 

0.0104 

Covariance 

Var (AGE = 5) 

oj 

0.0023 

0.0009 

Var (AGE = 10) 

O 2 ! 

0.0012 

0.0005 

Var (AGE = 15) 

O? 

0.0012 

0.0005 

Var (AGE = 20) 

°4 

0.0010 

0.0004 

Correlation 

P 

0.8664 

0.0573 


3-2). Except for the ARH (f) model (P = 0.101), no model 
was acceptable at the 5% level of significance compared 
with the unstructured covariance structure (P = 0.013- 
0.034). The AIC and BIC criterions also supported the 
ARH (1) covariance structure. Consequently, it was con¬ 
cluded that the model consisting of the quadratic mean 
structure and ARH (1) covariance structure best described 
the radial variation in wood density at breast height. Table 
4 shows the restricted maximum likelihood estimates and 
standard errors for all parameters in model 3-2 with ARH 
(1), denoted by model 4. Figure 3 shows the fitted average 
profiles of wood density for each spacing group. 

Effects of initial spacing on wood density 

In order to assess the effects of initial spacing on the radial 
variation in wood density, E-tests for the significance of sets 
of linear combinations of fixed effects in model 4 were 
applied. This could simplify the model by reducing the 
number of fixed effects. A first test assessed whether there 
was any quadratic age effect in the control or spacing treat¬ 
ments, and whether the slopes for the quadratic age effects 
differed between treatments. The null hypotheses tested are 
given by If: = /3 S = p 9 = 0 and H 2 : /5 7 = /) 8 = /$, (see Table 
5). Both hypotheses were rejected at the 5% level of signifi- 
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Fig. 3. Fitted average profiles for each spacing plot based on model 4. 
Thin solid lutes , plot A (300 trees/ha); broken lines , plot B (500 trees/ 
ha); thick solid lines , plot C (1000 trees/ha) 


Table 5. The null hypothesis and the P value resulting from the con¬ 
trast tests of fixed effects 



Null hypothesis 

p 

Hi 

A = A = ft = o 

0.003 

h 2 

ft = ft = ft 

0.035 

H 3 

ft = ft 

0.261 


ft = ft 

0.128 


ft = ft 

0.010 

h 4 

ft - ft - ft - o 

0.006 


ft - ft = ft 

0.050 

h 6 

ft = ft 

0.188 


ft = ft 

0.265 


ft = ft 

0.018 


cance (//,: P = 0.003; H 2 : P = 0.035). However, all pairwise 
differences (//,) were significant only for = /i, (P = 0.010). 
This indicates that the variation in wood density did not 
clearly correspond to the tree spacing intensity. Because the 
main purpose of this study was to investigate the effects of 
initial tree spacing, the mean structure with respect to the 
quadratic effect of age could not be simplified further. The 
same result was obtained for the tests of linear age effects, 
which were not zero and had slopes that differed between 
spacing treatments (H 4 : P = 0.006; H 5 : P = 0.050), although 
all pairwise differences ( H 6 ) were significant only for /f = 
P 6 (P = 0.018). These results suggest that wood density at 
breast height changes quadratically over time for each tree 
spacing treatment. 

The effects of initial tree spacing on wood density have 
been studied in a number of species. Kang et al. 22 reported 
in jack pine ( Pinus banksiana ) that wood density decreased 
with increasing initial spacing, and trees in the medium 
breast height diameter class were the most significantly 
affected. Johansson 20 found different results for Norway 
spruce {Picea abies ), whereby a decrease in initial spacing 
from 2.5 x 2.5 m (1600 trees/ha) to 1.5 x 1.5 m (4444 trees/ 
ha) had only a small effect, of little practical importance, on 
wood density. Larocque and Marshall 21 found similar results 
in red pine ( Pinus resinosa ) and suggested that closer 
spacing resulted in a faster increase in wood density from 


the pith outward. In the present study, wood density in trees 
in the control plot increased from 10 years of age, and the 
change was faster than that for trees in the treatment plots 
(Fig. 3). As Zobel and Sprague 26 defined, the juvenile wood 
zone is the area of rapid changes in properties near the pith 
and the mature wood is more uniform toward the bark. It 
is not possible to interpret the effects of initial spacing 
on the radial variation of wood density comprehensively 
through this short-term experiment; however, initial spacing 
may influence the age of transition from juvenile to mature 
wood. 

The growth and stem-form characteristics (Table 1) were 
compared among the various tree spacings. For all charac¬ 
teristics except total height (TH), highly significant differ¬ 
ences were found between the treatment groups (P = 
0.001-0.0014; data not shown). Larocque and Marshall 21 
stated that wood density increased more slowly with age, 
and the crown ratio decreased more slowly as tree spacing 
increased. A similar trend was found in this study (see Table 
1 and Fig. 3). The effects of initial tree spacing on the 
growth and stem-form characteristics were quite evident, 
but the effects on the wood density value were unclear 
despite the radial trends in wood density differing between 
the treatments. 

The model presented in this study is based on measure¬ 
ments at breast height and involved parameters including 
initial tree spacing and age effects. It is therefore necessary 
to extend this study to examine variation in wood density 
within the whole trunk. In addition, if it included parame¬ 
ters related to growth rate, the model will be useful for tree 
plantation managers in making decisions on when and how 
much to cut to achieve optimum timber quantity and 
quality. 
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